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Abstract 

Starting from the kinetic formulation of the hard thermal loop effective 
theory, we have (re) derived the collision terms for soft modes of order 
g 2 T log(l/<?) by averaging the statistical fluctuations in the plasma. 
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Over the past year, significant progress has been achieved in understanding the effective 
dynamics of soft excitations p ~ g 2 T in the non abelian quark-gluon plasma at high tem- 
perature. Bodeker has derived a diffusive field theory that involves only gauge fields with 
dynamics governed by a Langevin equation, in which the gaussian noise is only parametrised 
by the color conductivity, which shows a logarithmic dependence on g. As an application, 
the numerical value of the leading log behavior of the sphaleron rate, which measures baryon 
number violation in the hot standard model, has been computed ||. 

The crucial point in Bodeker's original derivation is the way in which the integration of 
hard ~ T and semi-hard ~ gT scales is done. The contribution of hard fields is easy to find. 
It is encoded in the local formulation of the hard-thermal-loop effective action J3],f|] leading 
to a closed set of collisionless kinetic equations and field equations for the semi-hard modes. 
The remaining integration of the semi-hard scales turns out to be a difficult task. Bodecker's 
main result is an effective kinetic equation with collision terms arising from the averaging 
over semi-hard fields. This equation has also been proposed by Arnold, Son and Yaffe || 
by analyzing the scattering processes between hard particles in the plasma. Very recently, 
a Boltzmann equation has been rigorously derived by Blaizot and Iancu |J, starting from 
the Dyson-Schwinger equations. For the QCD plasma, it coincides with the one previously 
obtained in Refs. fHH 

The main purpose of this paper is to give an alternative derivation d la Balescu-Lenard 
@ of the collision terms for the longwavelength color deviations of equilibrium. The starting 
point is the set of "microscopic" dynamical equations coming from the hard-thermal loop 
effective action describing the evolution of the collisionless plasma. In this way, collision 
terms arise as statistical averages of correlators of the plasma fluctuations at equilibrium 
H|9|]. A similar derivation starting from classical transport theory has been given recently 
by Litim and Manuel ||10|. It should be emphasized that the main point here, is to provide a 



derivation directly from the HTL effective action. As shown in Ref. p0| , a strictly classical 
derivation leads to a different value for the color conductivity because the value of the Debye 
mass computed from a Maxwell-Boltzmann distribution is different. 



f{x) 



9 J [2N c 5N(p,x) +N f 5n+(p,x) - N f 5n.(p,x)] . (2) 



Let us outline our calculations. Our construction relies on the local formulation of the 
HTL theory ||, which involves a set of coupled equations of motion for the mean fields 
and their induced current. The mean fields A a ^(x) satisfy the Yang-Mills equations with an 
induced current on the right hand side: 

[D„F*\ a = %, (1) 

where D u O(x) = d u O(x)+ig[A u , 0(x)\, A v = A a v t a and = [D^ D v ]/{ig). (The generators 
of the SU(N C ) gauge group in the fundamental representation are denoted by t a ; they satisfy 
[t a ,t 6 ] = if abc t c and tr(t a t b ) = S ab /2.) Furthermore, the induced color current j M is related 
to the color fluctuations in the gluon color densities 8N{p,x) = 5N a (p,x)t a and the quark 
color densities Sn±(p,x) = 5n±(p,x)t a by 

d 3 p 

In this equation, t> M = (l,v), v = p/p and Nf is the number of quark flavors. The system 
is closed by the non-abelian generalization of the Vlasov equations: 

[vD,5N(p,x)} -jvE^, (3) 

[vD,5n±(p,x)} =^v-E(x)^, (4) 

where N(p) = l/(exp(/3p) — 1), nip) = l/(exp(/3p) + 1) and E l a = F^°. The current is 
covariantly conserved. 

These equations contain the semi-hard degrees of freedom. In order to obtain a kinetic 
description for the soft modes, we decompose the fields into mean values and fluctuations: 

A^T a + 5K , (5) 
5N ^5N + 5N , (6) 
5n± — > 5n± + 5n± . (7) 

and average the previous microscopic equations. As a result, one obtains a set of equations 
for the average soft gauge fields A a ^(x), and for the mean values of the distribution functions 
IN, 5n±: 



D F 



v-D,SN(p,x) +gv-<E(x)> 
v ■ D, 5n±(p, x) ± gv • < E(x) > 



= + non-linear terms 
dN(p) 



dp 
dn(p) 
dp 



-ig v- < [5 A, 5N] > , 
—ig v- < [5 A, 5n±] > 



(8) 
(9) 

(10) 



where D u = d u + ig[A u , •] and < E J >= . Thus, the set of equations above includes 
the "collision integrals" which are determined by the correlators of the fluctuations. A 
noteworthy feature of the collision terms is their non-Abelian origin. A similar analysis 
starting with the abelian Vlasov-Maxwell equations for a QED plasma leads to no collision 
terms because of the linearity underlying this case. 

The main difficulty in computing the collision terms is caused by the nonlinearity of the 
equations for the fluctuations. However, the situation becomes significantly simpler for a 
plasma at high temperature with the coupling constant assumed to be small. In this case, 
we can neglect terms of higher order in the fluctuations and the following linearized 
equations for the fluctuations can be used: 



D v , 5F Uf * 
v ■ D, 6N(p, x 

v ■ D, 5n±(p, x) 



+ ig 

>, 

+ igv 
+ igv 



SA U ,F^ 



SA,SN 
SA, 5n± 



—g v • 5E(x) 
=F# v ■ 8E(x) 



dN(p) 

dp 
dn(p) 

dp 



(11) 
(12) 

(13) 



where 5F flu = D^bAy — DvbA^ and 5E l = 5F t0 . In order to obtain the explicit form of the 
kinetic equations we are looking for, it will suffice to linearize the collision terms in the mean 
values of the distribution functions, since the left-hand sides of Eqs. (p|)-(|IOD are already 
linear functions of them. Thus, the task is to compute the fluctuations up to first order in 
5N, 5n± and zero order in A. Moreover, I will only consider mean color distributions with 
no dependence in (t, x), which means that the results to be derived must be considered as 
zero-order terms in some derivative expansion. 

In order to implement the suggested programme, it is customary to use Laplace trans- 
forms, such as 
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5N a {u, k; p) = J d 3 x jf°° dt e- i( - k - x -^5N a {t, x; p) , (14) 
SA^u, k) = / d 3 x [°° dt e'^-^SAUt, x) , (15) 



In terms of these, when the Coulomb gauge is used, the first order correlator can be written 

as 

< M>,k)5iV 6 (u/,k'; P ) >«= l lf^E (P)^' < Ml(a;,k)^ m (a; / ,k / ) > 

Co> — K • p 

j^g2 jaran 



X 



A . 2 ^ , u , /v) < *4J»(", k)^^', k'; p) > (>' - 

/ (t^-Tp- K^V) + ^(PO " ^-(P)) , (16) 

with a similar expression for < 5A % a (u, k)8n±(uj', k'; p) >W. In this equation, e_i_(u;,A;) 
represents the transverse part of the dielectric function, defined by 

e ij {uj, k) = e\\(u, k)Uk 3 + e ± (u, k){5 ij - titi) 

Eqs. (jl6|) and (|17D deserve some comments: 

i) The prescriptions — > u;-M0 + , ct/ — > u/-M0 + are implicit in all denominators because 
the (double) Mellin's inversion formula for Laplace transforms requires that all singularities 
lie in the lower half u, to' planes. 

ii) the correlators over the initial conditions for t = 0, 



< 5N a (k; p)5N b (k'; p') > = (27r) 3 5(k + k')5 ab ((2tt) 3 5(p - p')f N (p) + fi N {k; p, p')) , (18) 

< &4(k; p)5n b ± (k'; p') > = (2vr) 3 5(k + k')5 ab ((2tt) 3 5(p - p')/ n (p) + /i„(k; p, p')) , (19) 

< <L4 l a (k)<L4f (k') > = (2vr) 3 <5(k + k')5 ab (8 ij - PP) fi A {k) , (20) 



have played an important role, since in Eqs. flTfiP and ( ]T7| ) only potential non-zero contribu- 
tions to the large time behavior have been written. This is so because, as proven in ||, when 
fN,n{p), fJ>N,n(k; p, p')) and /iyi(k) are smooth functions of the momenta, the only contribu- 
tions giving a non damped function of t and t' come from the functions /jv,n(p)- Thus, the 



initial conditions for 8 A, d8A/dt and the specific form of the two-particle correlation func- 
tions do not matter. All the information about the large time behavior is included in /jv,n(p) 
whose explicit form can be determined as follows. After substitution of Eqs. (|1~8|)-(|2"0|) into 
the appropriate correlators we find 



< 5A\ ( u, k)8Ai{u', k') >= (27r) 4 5(k + k')6{u + u')5 ab (S ij - titi) S(u, k) 

= (27r) 4 5(k + k')5{u + u')5 ab (5 ij - k r k j ) (21) 
g 2 (l-LU 2 /k 2 ) 1 



J o °° dpp 2 (AN 2 J N (p) + 2N 2 f n (p)) . (22) 



I — Lu 2 e±(u, k) + k 2 \ 2 Ank Jo 
Then, a comparison with the prediction of the fluctuation-dissipation theorem relating the 
correlator of field fluctuations with the imaginary part of the retarded propagator, 

S , (a;,A;) = (l-e- /Sw )- 1 2Im— — — - , (23) 

V ; V ; k 2 -(LU + iO + ) 2 e ± {LU + iQ+,k)' V ; 

reveals the right values for the initial correlators, 

f N (p) = N' 1 N(p)(l + N(p)) = —N~ l TN'(p) , (24) 
f n (p) = 2Nj 1 n{p){\ - n{p)) = -2Nj x Tn'{j>) , (25) 

where the factors iV" 1 and 2NJ 1 compensate similar factors in the color current. These 
results for the initial correlators could also be established from the general properties of the 
fluctuations of an ideal gas [XT]. Up to the color and flavor factors, Eqs. ( p4|) and ( p5|) are 
exactly the equations for the mean square of the fluctuations presented in Ref. ||ll|f|. 

Now, following the Ref. @, other correlators can be easily computed. For example, at 
leading order in g we have 

< 5Ai(u, k)5N b (u', k'; p) > = (2?r) 4 5(k + k')5(u + J) Ang 5 ab - k l k j ) p> 

x ^) ( i + ^)) (26) 

k 2 — u 2 e±(uj, k) 

< 5A l a {u, k)5n b ± {u', k'; p) > = ±(2vr) 4 5(k + k')8(u + uf) 2ixg 5 ab (5 ij - Pk j ) p> 

n(p)(l — n{p)) 



k 2 — u 2 e±(u, k) 



6{u-k-p). (27) 



1 Specifically, Eqs. (113,3) and (113,4) of Ref. Q. 
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iii) The dominant scattering processes in the collisions correspond to i-channel exchange 
of quasistatic magnetic gauge bosons, which means that the longitudinal part of the interac- 
tion can be ignored. Also, as the relevant regime to be considered in transverse propagators 
is to <C k, we can make the following replacements 

37riu; 2 



e±(u,k) 



4uk ' 
S(uj,k)~^-8(uj) 



k 2 p 
2 



5(u) 



(28) 
(29) 
(30) 



(-k 2 + u 2 e ± (uj,k)) 2 3ku 2 
with the plasma frequency to 2 = g 2 {2N c + Nf)T 2 / 18. 

Finally, the collision terms easily emerge from all of this by integration over (u, k). Their 
explicit form turns out to be 



-igv- < [5 A, SN] a >= -j g 

d 3 p' (p • p') 2 



4^2 



X 



-igv - < [SA, 5n±] a >= -j g 



A2 



2N c SN a (x, p') + N f 8n + {x, p') - N f 5n_(x, p' 



5n±(x,p) ± ^jl n \P) 



x 



d 3 p' 



IP P 



(27T) 



r\2 

J== (2N c SN a {x, p') + N f Sn + {x, p') - N f 5n a _{x, p' 



:p • pO 



with the damping rate 7 5 for the hard thermal gauge bosons 

g 2 N c T 



la 



Air 



log(l/<?). 



If we combine the different color distributions functions in the form |3||| 
1 f dp Amp 2 / —— a 



W a {x,p) 
we obtain 



3u 2 



(27T) 



(2A^ c 5AT a (x, p) + N f &n + (x, p) - N f &^_(x, p] 



(31) 



(32) 



(33) 



(34) 



u • D, W(a;, p)] - v ■ < E a (x) >= -SC[W} a , 

5C[I^] a = 7 [ p) - - / dO, (P • P S = v^°( X) p/) 

•/ -a- (p. p') 2 



(35) 
(36) 



from which the color conductivity o = oj'i/^g follows [|1|,|5|,|12 



In conclusion, I have derived from the microscopic HTL equations the explicit form of 
the collision terms recently proposed by Arnold, Son and Yaffe || and Blaizot and Iancu || 
which are relevant for Bodeker's effective theory As the authors of the Ref. [|T0]] have pointed 
out in a strictly classical context, this approach is simple. It is based on the properties of the 
plasma fluctuations at equilibrium. In this sense, this approach would be closely connected 
with linear response theory based on "Kubo" formulas expressing the transport coefficients. 
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